S-ROCK: Chebyshev Methods for Stiff Stochastic Differential Equations

We present and analyze a new class of numerical methods for the solution of stiff stochastic differential equations (SDEs). These methods, called S-ROCK (for stochastic orthogonal Runge–Kutta Chebyshev), are explicit and of strong order 1 and possess large stability domains in the mean-square sense. For mean-square stable stiff SDEs, they are much more efficient than the standard explicit methods proposed so far for stochastic problems and give significant speed improvement. The explicitness of the S-ROCK methods allows one to handle large systems without linear algebra problems usually encountered with implicit methods. Numerical results and comparisons with existing methods are reported.